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ABSTRACT 

In this paper the plane elasticity problem for two bonded half planes containing a crack 
perpendicular to the interface is considered. The primary objective of the paper is to study the 
effect of very steep variations in the material properties near the diffusion plane on the 
singular behavior of the stresses and stress intensity factors. The two materials are, thus, 
assumed to have the shear moduli /z 0 and /i Q exp(/?x), x=0 being the diffusion plane. Of 
particular interest is the examination of the nature of stress singularity near a crack tip 
terminating at the interface where the shear modulus has a discontinuous derivative. The 
results show that, unlike the crack problem in piecewise homogeneous materials for which the 
singularity is of the form r~ a , 0<a<l, in this problem the stresses have a standard square-root 
singularity regardless of the location of the crack tip. The non homogeneity constant /? has, 
however, considerable influence on the stress intensity factors. 

1. Introduction 

From the viewpoint of mechanical failure in diffusion bonded materials and other 
adhesively bonded structures one of the most important problems is, along with the fracture 
characterization of the interfacial zones, a proper mechanical modelling of these zones and the 
development of suitable techniques for solving the related crack problems. The significance of 
the problem lies in the fact that (i) regardless of the mechanism of binding at the atomic level, 
in many cases there is always a thermodynamically stable and readily distinguishable region 
between the two homogeneous materials, (ii) the interfacial regions are generally locations of 
higher concentration of stresses and microflaws, and (iii) in most material pairs the fracture 
toughness and the subcritical crack growth resistance of the interfacial zone tend to be lower 
than that of adjacent homogeneous materials. The standard continuum solid mechanics model 
for bonded materials has been the usual piecewise homogeneous medium with an ideal interface 
across which the physical properties of the medium exhibit a jump discontinuity. For many 
applications this model is quite adequate. However, there are two groups of fracture related 
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problems in which the model leads to some physically anomalous results and may require 
modifications. The first is the complex power singularty at the tip of an interface crack giving 
oscilations in stresses and crack opening displacements very near the crack tip (e.g., 
Muskhelishvili (1953a), Erdogan (1965), England (1965), Rice and Sih (1965)). A crack tip 
contact model has been proposed to eliminate these oscillations within the context of the linear 
elastic theory, Comninou (1977), Atkinson (1977). 

The second problem deals with a crack intersecting the interface. In- this case it is 
known that the stress state has a power singularity of the form r‘ a with a^l/2 and 0<a<l, 
Cook and Erdogan (1972), Erdogan and Cook (1974), Erdogan and Biricikoglu (1973). The 
physically objectionable feature of this result is that because of the non square- root singularity 
the propagating crack would not remain self-similar, and consequently the standard energy 
balance based fracture mechanics techniques would not be applicable to problems in bonded 
materials that involve cracks intersecting the interfaces. 

The ideal interface model may be improved by considering the interdiffusion nature of 
the bonding process. In some cases the resulting reaction zone may be very thin and in others 
it may consist of intermediate phases or intermetallics with mechanical properties substantially 
different than the adjacent materials. In either case the microprobe line scans seem to indicate 
that in the interfacial region the percent distribution of atoms from one phase to the next is 
very steep but, nevertheless, continuous. Thus, in diffusion bonded materials one could replace 
the ideal interface requiring a jump discontinuity in physical properties across the boundary by 
a thin nonhomogeneous layer with highly steep variations in the mechanical properties. One 
may further assume that the elastic parameters of the composite medium are continuous 
functions of the space coordinates throughout the medium having discontinuous derivatives 
along the boundaries of the interfacial region. The analytical questions raised by this material 
model concern the effect of the material non homogeneity, particulary of the “kinks” in the 
distribution of the elastic parameters on the nature of the crack tip singularities. The crack 
problem for an infinite nonhomogeneous medium was considered by Delale and Erdogan (1983) 
who showed that as expected the standard square-root singularity at the crack tips is 
maintained if the material parameters are continuous and differentiable functions. 

Regarding the physical anomalies that follow from the piecewise homogeneous material 
model, the first question deals with the nature of singularities for cracks lying along the kink 
lines of the elastic property distribution. This question was recently considered by Delale and 
Erdogan (1988a) who showed that for these “interface” cracks not only the square- root 
behavior of singularity but also the angular distribution of stresses very near the crack tips 
turn out to be identical to that found for the homogeneous materials. Further results for a 
crack in a nonhomogeneous interfacial zone are given by Delale and Erdogan (1988b). The 
second basic question in introducing the nonhomogeneous interfacial zone model is concerned 
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with a crack intersecting the interface, specifically with the nature of stress singularity for a 
crack terminating at the kink line of the elastic property distribution. The problem under 
anti- plane shear loading was considered by Erdogan (1985) who showed that the square-root 
nature of stress singularity remains unaffected by the material nonhomogeneity provided the 
shear modulus is a continuous but not necessarily a differentiable function near and at the 
crack tip. The behavior of stress singularity in the corresponding plane strain problem is one 
of the remaining unresolved questions that will be considered in this paper. 


2. Formulation of the Crack Problem 

Consider the plane elasticity problem shown in Fig. 1 where a homogeneous half plane 
containing a crack perpendicular to the boundary is bonded to a nonhomogeneous half plane. 
For the purpose of examining the nature of crack tip singularity and for studying the effect of 
material nonhomogeneity on the stress intensity factors, it will be assumed that the elastic 
properties of the two materials are given by 


/ \ /?x 

P l(x) = /*oe 


, p 2 = Ho , K X = K 2 = K 


(i) 


where fi is the shear modulus, /c=3-4i/ for plane strain and /c=(3-t/)/(l-f v) for generalized 
plane stress, v being the Poisson’s ratio. With (1) the Navier’s equations may be expressed as 


V 2 n . , _2_ a , iZh ,0/9^4. = n 

i k- 1 0x 2 dxdy' ^ dx ’ 


dy- 


V 2 % 


2 /d 2y i 
+ — 1 (^2 + 


Sp + ^ + £) = ° ’ 


(2a, b) 


where Uj and v j? (i=l,2) are the x and y components of the displacement and /3 ^ 0 for i=l or 
x<0 and /?=0 for i=2 or x>0 (Fig. 1). Referring to Fig. 1, it will be assumed that y=0 plane 
is a plane of symmetry for the external loads as well as the geometry. It will also be assumed 
that the problem in the absence of the crack under the actual loading conditions has been 
solved and ^2y y (x,0) has been determined. Thus, the important results concerning the crack 
problem may be obtained by considering only the perturbation problem for y>0 in which the 
normal traction p(x)=-<72 yy (x>0) is the only external load. Equation (2) must then be solved 
under the following conditions: 


Uj(0,y) = u 2 (0,y) , v 1 (0,y) = v 2 (0,y) , 0<y<oo , (3a, b) 

<Tixx(0,y) = <7 2 xx(°>y) , <^ixy(0,y) = <7 2 xy(°’y) » 0<y<oo , (4a, b) 
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v i(x,0) = 0 , <7 lxy (x,0) = 0 , -oo<x<0 , 


(5a, b) 


<T 2xy( x >®) — 0 ) 0<x<oo , 


( 6 ) 


v 2 (x,0) = 0 , 0<x<a , b<x<oo , 


°2yy( x »+°) = P( x ) > a<x<b . 


(7a, b) 


The solution of (2) may now be written as 

f °0 

ii(x,y) = | f(x,a)cosayda , 

J o 

/■OO 

v l( x >y) = | g(x,a)sinayda , -oo<x <0 , 
j o 

r°° r°° 

u 2 ( x »y) = I f 1 (x,a)cosayda + ± f 2 (y,a)e‘ iax da , 

Jo w J-« 

o r°° r°° 

v 2 ( x »y) = I g 1 (x,a)sinayda + *L g 2 (y,a)e' iax da , 0 <x<oo , 0 <y<oo.( 9 a,b) 

j O J-oo 


(8a, b) 


By substituting from (8) into (2) it can be shown that 

A n : x 4 n=x 


f(x,a) = £ mjDj(ar)e J , g(x,a) = £ Dj(a)e J , (10a, b) 

j=l j=l 

where D 1 ,...,D 4 are unknown, n lv ..,n 4 are the roots of the characteristic equation 


n 4 + 2/?n 3 + (/? 2 -2a 2 )n 2 - 2/?a 2 n + (^)/?V + a 4 = 0 , 


( 11 ) 


and m lv ..,m 4 are given by 

nj 2 +/?nj-(l+ao)a 2 


m j = 


(a c n j +/?)a ’ 80 
By rewriting (11) in the following form 


_ 2 _ j — 1 4 

k -1 ’ J ~ lv ”’ 4 


( 12 ) 


(n^+/?n-a 2 ) 2 + (^)/? 2 a 2 = 0 , 




the roots may be obtained as 
1 


n l = ^(-/?+Rcos 7 ) + ^ Rsin 7 , n 2 = fi 1 , 


( 13 ) 


(14a, b) 
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(15a, b) 


n 3 = |(-/?-Rcos 7) + ^ Rsin7 , n 4 = n 3 , 


/ 4/?a v3T 

R = [(/> 2 +4« 2 ) 2 + 16^a 2 (^)] 1/4 , t.„2 T = 


(16a, b) 


where we note that R is always positive and 


Re(n 1 ,n 2 ) > 0 , Re(n 3 ,n 4 ) < 0 . 


(17) 


obtain 


Observing that for the material 1 x<0, by using the regularity conditions at infinity we 

pOO 2 

u 1 (x,y) = | I £ mj(o)Dj(a)e j cosayda , 

Jo 1 

v x (x,y) = f S Dj( a ) e ^ ^ sinayda , -oo<x<0 , 0<y<oo , (18a, b) 

Jo 1 


r lxx 


lyy 


(x,y) = ^T” * I ^ [(l+")mjnj + (3-K)a]Dje njX cosayda 


/?x -oo 2 n-x 

(x,y) = £^5- | £ [(l+*)a + (3-«)mjnj]Dje J cosayda , 

Jo 1 

/? f 00 2 n x 

(x,y) = p 0 e /?X | £ (nj-amj)Dje J sinayda , -oo<x<0 , 0<y<oo , (19a-c) 

Jo 1 


Similarly, substituting from (9) into (2) and letting /?— 0, we find 

u 2 (z,y) = I [ [A x +(x + a)A 2 ]e* ax cosayda 
Jo 

f OO 

v 2 (x,y) = | (A 1 +A 2 x)e' ax sinay 
Jo 

• i f” (j£j + y) , 0 <x<oo , 0 <y< 


(20a, b) 
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2 ]r o <7 2 xx( x >y) = - I | i aA i + ( ax + ^-)A 2 ]e' ax cosayda 
i f°° A<*\ \* -laly-iax 

+ i (W - «y)A 3 e da » 

J -oo 

fOO 

2 ~ <T 2yy( x ’y) = I j ( aA l + ( ax + ^ 3 )A 2 ]e- ax cosayda 

: f 00 .lal . . -laly-iax, 

+ ^ (W + ay)A 3 e da , 

J -00 

f OO 

2 ~ ^ 2 xy( x >y) = - 1 j t aA l + ( ax + ^) A 2 ] e ' aXsin£ *yda 

- ^ |a|yA 3 e * aX da , 0<x<cx> , 0<y<oo , (21a-c) 

J o 

where A 1 , A 2 , A 3 are unknown functions. 

It may be seen that because of the assumed symmetry in (8) and (9), the conditions (5) 
and (6) are identically satisfied. The remaining five conditions (3), (4) and (7) may then be 
used to determine the unknowns D x , D 2 , A x , A 2 , and A 3 . 


3. Derivation of the Integral Equation 

To reduce the problem to an integral equation we first introduce 

4>(x) = v 2 (x,0) , a<x<b . (22) 

From (22) and (20b) it then follows that 



The four homogeneous conditions (3) and (4) may now be used to express D 1? D 2 , and A 2 
in terms of <j>, which may then be determined by using (7b). The expressions of D x , D 2 , A 1 
and A 2 are given in Appendix A. 

In considering the boundary condition (7a), from (21b) and (23) the contribution of the 
second term on the right hand side of (21b) may easily be evaluated as follows: 


f (!gl + oy )A 3 (c,)e-l“'^ i “ X d ( , 

j -oo 


—/ — am 
n«+l) y->+0 




f+y 2 
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- it o f *(*)<** (oa\ 

”*(*+l)J a *- x * . . ( j 

Evaluating the contribution coming from the first term on the right hand side of (21b) by 
making use of the results given in Appendix A, the condition (7b) may now be expressed as 


if [ + K(x,t)]*(t)dt = P(x) , a<x<b , 

* fi. 


(25) 


where the kernel K is given by 

-a(t+x) 


K(x,t) = f h(a,x,t)e ” v * ' "'da , 

Jo 

h(a,x,t) = h^a) + xh 2 (a) + th 3 (a) + txh 4 (a) , 


(26) 

(27) 


hi(a) = (*c+l)[ap 3 (a) + ^ p 4 (a)] , 


h 2 (a) = («+l)ap 4 (a) , 

M°0 = (<c+l)[aq 3 (a) + ^ q 4 (<*)] , 

h 4 (a) = (K+l)aq 4 (a) • (28a-d) 


The functions p 3 , p 4 , q 3 and q 4 may also be found in Appendix A. 

Before proceeding with the examination of the asymptotic behavior of the kernel K(x,t), 
two special cases may be considered. The first is the case of /?= 0 for which the problem is one 
of a homogeneous infinite plane containing a single crack. For /?— ► 0 it may easily be shown 
that A— m 2 , B— ►O, C— *-l, D— *0 and, consequently K(x,t) vanishes giving (see Appendix A) 

i j a i^r dt = p ( x ) » a < x < b ( 29 ) 

which is the expected result. In the second special case we let /?=oo which, from Fig. 1 and 
p x (x) = p 0 exp(/?x), (x<0), is seen to correspond to p 1 =0. The problem would then be one of 
a homogeneous half plane having a crack perpendicular to the boundary. From (12), (14), 
(16), (A15) and (A16) it can be shown that for /?— *-oo 


A — 0 , B 


a. 


3-k 

1+K 



( 30 ) 


Thus, by substituting from Appendix A and (30) into (28) and (27) we obtain 
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h(a,t,x) — ► (2-ax-3at+2a 2 tx) . 
From (26) and (31) it then follows that 


( 31 ) 


K(x,t) = - 4* + 


6x 


Azl 


(t+x)-* (t+x) ; 


(32) 


which is the expected half plane result (Cook and Erdogan, 1972). 

In piecewise homogeneous materials having the same crack geometry the kernel K has a 
form quite similar to (32), namely (Erdogan and Biricikoglu, 1973) 


K (x t) = + b2 * .. , Ja*. 

( ’ ; t+ x + (t+x) 2 + (t+x) 3 


(33) 


where b l9 b 2 , and b 3 are material property dependent constants. Note that if the crack tip 
touches the interface, i.e., for a=0, K would become unbounded as x and t approach zero and 
hence would contribute to the singular behavior of the solution. It is then clear that in order 
to examine the singular nature of the solution of (25), the asymptotic behavior of K must be 
investigated. Since the function h in (26) is integrable at a=0 and bounded elsewhere, any 
possible divergence in K must be due to the asymptotic behavior of h as a— *oo. Referring now 
to Appendix B, for a— >oo we obtain 

h(a,x,t) -► hoo(a,x,t) = + (d 2 + ^)t + (d 4 a + d 5 + ^)tx , (34) 


where 


<*1 



d 2 — -d 4 


£ 

1 + K 


, d 3 — -d 5 


( l -«)/? 2 

2 ( l +/ c ) 2 




/ c ( 3 -/ c )/? 3 
4 ( l +/ c ) 3 ' 


(35) 


By subtracting and adding hoo from and to h in (26) and by evaluating the integrals, it 
may easily be shown that 


k(x,t) = k 1 (x,t) + k 2 (x,t) + k 3 (x,t) , 


(36) 


k 1 (x,t) = f (d 2 t+d 5 tx+d 4 atx)e a ^ t_ *" X ^da 

Jo 


= d 


2 t+x 


+ d c + 

+ t+X + 


tx 


( t + x )" 


(37) 
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k 2 (x,t) = -(d 1 +d 3 t+d 6 tx)Ei[-a 0 (t+x)] , 


(•“O 

k 3 (x,t) = [h(a,x,t) - (d 2 t+d 5 tx+d 4 txa)]« 

Jo 

+ [ [h(a,x,t) - h«(a,x,t)]e “^"^da , 

Ja n 


-a(t+x) 


da 


( 38 ) 


where a Q is any positive constant 2 and the exponential integral is defined by 
Ei[-a 0 (t+x)] = -[ ^ e a ^ +X ^da 

Ja 0 

= C 0 + loga Q + log(t+x) + £ ^ Q ° , 

k=l k ^ k -J 

a 0 > 0 , t+x > 0 , (39) 


and C 0 = 0.5772156649 is the Euler’s constant. 

By substituting from (36)-(38) into (25) we obtain 

J l Jtfe + ( d 4' d 2) tfe ' d 5 - d 4 ^j2^ (t ^ dt 

+ \ J [d 2 +d 5 x+k 2 (x,t)+k 3 (x,t)]<^(t)dt = p(x) , (a<x<b) . (40) 

For a>0 (40) is an ordinary singular integral equation the dominant part of which has a 
simple Cauchy singularity. Its solution may, therefore, be expressed as (Muskhelishvili, 1953b) 


*(t) = 


F(t) 

^ (t-a)(b-t) 


, a<t<b , 


(41) 


where F is a bounded unknown function. For a=0 the kernel of the second term in the left- 
hand side of (40) can still be treated as a Fredholm kernel and would not have any 
contribution to the singular behavior of <£(t). To examine the first term we let 


<t>{ t) = — V —f\ , 0<t<b , 0<Re(u>,r/)<l , 
t w (b-t) ' 

where F again is a bounded function. Defining the sectionally holomorphic function 


2 In numerical analysis rapid convergence is obtained by selecting 0<a o <l. 
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( 42 ) 



it can be shown that (Muskhelishvili, 1953b) 


V>(z) 


F(0)e i7ru) ± _ F(b) i 
b^sinxw z w b w sin5r»7 (z-b)^ 


+ V’o(z) , 


(44) 


where V’o(z) may have at most weaker singularities than z‘ w and (z-b) ** at z=0 and z=b, 
respectively. From (43) and (44) it can be shown that 


1 r b <000 J, _ F(0)cot7TW 1 
Mo ^ b’ 7 x w ' 


F(b)cotirr; 

b° 


^~fj + V’i(x) , 0<x<b 


(b-x) 


(45) 


where at the end points the behavior of is the same as that of o- From (43) and (44) it 
can also be shown that near x=0 we have 


if dt = V-(-x) = + V> 2 (x) , 0<x<b 

Jo ~ b sinxw x 


(46) 


where at x=0 tfi 2 again may have at most a weaker singularity than x*^. From (46) it follows 
that 


1 

7T 


1 

7T 


f b x k (ft(t) 

Jo t+x 


f b X 2 ^(t) 

Jo (t+x ) 2 


dt 


dt 


X k ‘ w + x k V> 2 (x) , k=l,2,..., 0<x<b , 

b sin iruf 


^F(O) i-v 
u r/ • 

b sinTTUJ 


X 2 ^ V’2(x) , 0<x<b . 


(47) 


(48) 


If we now substitute from (45), (47) and (48) into (40), first multiply both sides by 
(b-x) 7 * and let x=b then by x u and let x=0, and observe that F(0)^0, F(b)540, we find 

cot iru> = 0 , cot 7rjy = 0 . (49a, b) 


giving w = T) = 1/2. This is an important result indicating that the solution of the problem 



has a square-root singularity regardless of the location of the crack tip and the solution has the 
form as given by (41) for a=0 as well as a>0 3 . 


4. The Solution and Results 

To solve the integral equation (25) we make a change in variable and define the 
following normalized quantities: 




(50) 


= m 

] (t-a)(b-t) ] 1-r 2 


a<t<b , -l<r<l 


(^)K(x,t) = k(s,r) , 


(51) 

(52) 


^ P( x ) = Q( s ) 


(53) 


Equation (25) then becomes 

I j 1 [A + k(r,s)] dr = Q(s) , -1 <s<1 . (54) 

From (7a), (22) and (23) it is clear that the integral equation must be solved under the 
following single-valuedness condition: 

f b 

<£(t)dt = 0 (55) 

J a 



(56) 


3 Based on this and the previous results found by Delale and Erdogan (1983,1988a), one 
may now conjecture that in plane elasticity problems the material nonhomogeneity would have 
no influence on the nature of crack tip singularity regardless of the location and the orientation 
of the crack as long as the material parameters are continuous, but not necessarily 
differentiable, functions of the coordinates. This means that for a crack intersecting the 
interfacial region at an arbitrary angle, the complex power singularity predicted by the 
piecewise homogeneous material model for certain materials and crack orientations would also 
disappear and the square-root singularity at the crack tips would be the only singularity the 
solution may have. 
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Equation (54) is solved by substituting 


G(r) = £ a„T n (r) ( 57 ) 

o ' 

and by using a simple collocation scheme to reduce the resulting functional equation to an 
infinite system. First, from (56) and the orthogonality condition of Chebyshev polynomials T n 
it may be seen that 


ao = 0 . 


Thus, from (54), (57) and (58) by truncating the series and by using 


1 f 1 Tn(r)dr _ J 

J-l (r-s)Vl-r 2 \ 

N , rl 

£ a n[U n . 1 (sj) + i 
1 J -1 


0 , n = 0 


Up_j(s) , n — 1 , 2 ,..., -1<s<1 , 


Tn(r)k(r,Sj) 
■> 1 -r 2 


dr] = Q(sj) , j=l,. 


■»N 


(58) 


(59) 


(60) 


where U n is the Chebyshev polynomial of the second kind. The collocation points Sj may be 
selected as follows: 


Tn(sj) = 0 , Sj = C os[(2j-l)^] , j=l,...,N . 


(61) 


Since (25) gives <r 2y y(x,0) outside as well as inside the crack, substituting from (41) into 
(25) a simple asymptotic analysis would show that near the crack tips we have (see, for 
example, Erdogan, 1985) 


<r 2 y 2 (x, 0 ) 


iifo F(b) 

(b-a)(x-b) 


+ ff l(x) . 


<r 2yz( x, °) 


4po 

l+/c 


F(a) 

| (b-a)(a-x) 


+ ^ 2 ( x ) 


(62a, b) 


where u ^ and <r 2 represent the part of the solution that is bounded at x=b and x=a, 
respectively. Thus, from (51) and (62) the mode I stress intensity factors defined by 
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k l( a ) = j 2(a-x) <r lyy (x,0) , k^b) = lim b \ 2(x-b) £r lyy (x,0) 


(63) 


may be expressed as 4 


M») = G(-l)^ , k l( b) = ■ G(l)^ 


(64a, b) 


The main results obtained in this study regarding the stress intensity factors are shown 
in Table 1, where it is assumed that 


= c , = d > p(x) = -<r 


(65) 


The uniform crack surface pressure a used in the example is realistic if the composite medium 
is subjected to a constant displacement loading in y direction away from the crack region. In 
the examples considered it is further assumed that v = 0.3 and the plane strain conditions 
prevail. By examining the results given in the table a number of observations may be made. 
First from Fig. 1 and /^(x) = /i Q ex p(/?x) we note that for /?>0 medium 2 is stiffer than 
medium 1, hence the stress intensity factors would be expected and are seen to be greater than 
(r-fc, the value for the corresponding homogeneous material having the same crack. For /?<0 
material 1 is the stiffer medium and consequently the stress intensity factors, as expected, turn 
out to be less than <r^ c. These trends are also observed for varying /?, that is for any crack 
geometry as j3 increases, the stress intensity factors increase. For f3— mx> the results shown in 
the table agree with that of a half plane containing a uniformly pressurized crack of length 2c 
perpendicular to the boundary. As the distance d increases, the stress intensity factors tend to 
tf'fc, the homogeneous plane solution. 

Figures 2-4 show the variation of the normalized stress intensity factors defined by 


k(.) = bW , Kb) = bW 


( 66 ) 


with the crack distance d and for various values of the nonhomogeneity parameter /?. For the 

4 The expression (64a) is valid also for the limiting case of a=0. In fact (3b) and (4a) 
implies that du/dx is continuous across x=0. Since <9v/<9y, /i and k are also continuous, this 
means that ^ lyy (-0,y) = ^ 2yy (+0,y), y>0, implying that the stress singularity will be 
controlled by that of </>(x) and the values of /i and k at x=0, y=0. 
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specific crack geometry a=0, the effect of /? on the normalized stress intensity factors are 
shown in Fig. 5. In this figure the absissa is selected in such a way that the results include the 
full range of /?, -oo</?<oo. As seen from the figure for ► oo the stiffness of medium 1 
approaches zero, k(a)-+oo and k(b) goes to the edge crack solution, 1.5861 marked on the 
figure. Similarly, for /?— »oo medium 1 becomes a rigid half space, k(b) goes to the appropriate 
limit shown in the figure and k(a) tends to zero. 

Finally, Fig. 6 shows some examples giving the normalized crack surface displacement 
for the limiting case of a=0 defined by 


V(x) 


v 2 (x,+0) 

(l + K)cr C /4/i 0 


(67) 


The dashed line shown corresponds to fi=0 or the homogeneous medium. Again, as expected 
the displacement is skewed toward or away from the interface depending on the sign of f3 or 
the relative stiffnesses of materials 1 and 2. 
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Table 1. Normalized mode I stress intensity factors 
for a crack of length 2c subjected to constant crack 
surface pressure a f I/-0.3, the case of plane strain. 


pc 

V a) 

k x (b) 

k^a) 

k x (b) 

k L (a) 

k x (b) 

aj" c 

crj~c 

aj c 

<7,fc 

ajc 

aj"c 


d/c-l . 0 

d/c-l . 1 

d/c-l. 25 

0.1 

1.0172 

1.0087 

1.0138 

1.0079 

1.0112 

1.0070 

0.5 

1.0820 

1.0377 

1.0638 

1.0337 

1.0499 

1.0291 

1.0 

1.1546 

1.0653 

1.1163 

1.0572 

1.0880 

1.0482 

2.0 

1.2791 

1.1042 

1.1980 

1.0884 

1.1426 

1.0721 

5.0 

1.5586 

1.1665 

1.3472 

1.1334 

1.2286 

1.1033 

10.0 

1.8785 

1.2138 

1.4692 

1.1620 

1.2870 

1.1211 

20.0 

2.3105 

1.2562 

1.5745 

1.1825 

1.3298 

1.1327 

50.0 

3.103 

1.3029 

1.6699 

1.1983 

1.3625 

1 . 1407 

-►CO 

00 

1.5861 

1.7585 

1.2108 

1.3875 

1 . 1464 


d/c-l . 5 

d/c-2 . 0 

d/c-5 . 0 

0.1 

1.0086 

1.0059 

1.0059 

1 . 0044 

1.0019 

1.0017 

0.5 

1.0367 

1.0235 

1.0234 

1.0167 

1.0057 

1 . 0049 

1.0 

1.0621 

1.0378 

1.0372 

1.0255 

1.0076 

1.0064 

2.0 

1.0950 

1.0543 

1.0530 

1.0347 

1.0091 

1.0076 

5.0 

1.1399 

1.0738 

1.0711 

1 . 0443 

1.0103 

1.0085 

10.0 

1.1662 

1.0839 

1.0802 

1.0488 

1.0108 

1.0088 

20.0 

1.1833 

1.0900 

1.0855 

1.0513 

1.0110 

1.0090 

50.0 

1.1952 

1.0940 

1.0890 

1.0528 

1.0111 

1.0091 

-KC 

1.2035 

1.0967 

1.0913 

1.0539 

1.0112 

1.0092 



d/c-l . 0 

d/c-l . 1 

d/c-l. 25 

-0.1 

0.9835 

0.9918 

0.9867 

0.9925 

0.9893 

0.9933 

-0.5 

0.9291 

0.9676 

0 . 9440 

0.9707 

0.9557 

0.9745 

-1.0 

0.8778 

0.9484 

0.9052 

0.9537 

0.9264 

0.9601 

- 2.0 

0.8057 

0.9264 

0.8529 

0.9345 

0.8888 

0.9443 

-5.0 

0.6880 

0.9008 

0.7750 

0.9127 

0.8374 

0.9269 

-10.0 

0.5963 

0.8879 

0.7228 

0.9019 

0.8066 

0.9184 

-20.0 

0.5122 

0.8801 

0.6840 

0.8955 

0.7860 

0.9135 

-50.0 

0.4181 

0.8749 

0.6534 

0.8913 

0.7712 

0.9103 

— ► - 00 

0. 

0.8710 

0.6283 

0.8883 

0.7601 

0.9081 


d/c- 

■1.5 

d/c- 

■2.0 

d/c- 

■5.0 

-0.1 

0.9918 

0.9944 

0.9944 

0.9958 

0.9982 

0.9984 

-0.5 

0.9671 

0.9791 

0.9789 

0.9850 

0.9949 

0.9956 

-1.0 

0.9469 

0.9680 

0.9675 

0.9779 

0.9933 

0.9944 

-2.0 

0.9227 

0.9563 

0.9553 

0.9710 

0.9921 

0.9934 

-5.0 

0.8929 

0 . 9440 

0.9422 

0.9642 

0.9912 

0.9927 

-10.0 

0.8768 

0.9382 

0.9360 

0.9613 

0.9908 

0.9924 

-20.0 

0.8669 

0.9349 

0.9324 

0.9596 

0.9906 

0.9923 

-50.0 

0.8601 

0.9328 

0.9300 

0.9585 

0.9905 

0.9922 

— ► - OO 

0.8553 

0.9313 

0.9284 

0.9578 

0.9904 

0.9922 
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h( a) 


C/So=. 1) 


( /So-. 5 ) 


(a) k(b) 


The variation of the normalized stress intensity factors with the crack distance, 
d=(b+a)/2, c=(b-a)/2, k(a)=k 1 (a)/ < r^c, k(b)=k 1 (b)/<rfa, the case of plane strain, 

i/=0.3, /?c=0.1 and /?c=0.5. 
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(00=5.0) 


(a) /c(b ) (0c=-5.O) 


I . I I I I L 


N 


(Bo=10. ) 


(a) (0o=-1O.) 


The variation of the normalized stress intensity factors with the crack distance, 
d=(b+a)/2, c=(b-a)/2, k(a)=k 1 )/ < rf£, k^^b )/aU, the case of plane strain 

i/=0.3, /?c=5 and /?c=10. 






1 /pc pc 


Fig. 5 The variation of the normalized 
constant /? for the special crack 


2 ; 



The normalized crack surface displacement for the special crack geometry a— 0 and 
for various values of /?, V(x)=v 2 (x,+0)/[(1+k)(o , c/4/i o )]. 




APPENDIX A 


The expressions of the functions D 1? D 2 , Aj and A 2 

By substituting from (18), (19b, c), (20), (21b, c) and (23) into (3) and (4), after rather 
straightforward algebra we obtain 

D x (a) = [ [ci(a,t) + ic 2 (a,t)]e' at ^(t)dt , 

J a 

D 2 (ct) = [ [ci(a»t) - ic 2 (ar,t)]e' at ^(t)dt , 

J a 

A x (a) = [ c 3 (a,t)e’ a V(t)dt , 

J a 

A 2 (a) = f c 4 (a,t)e* a V(t)dt , (Al-4) 

J a 

where 

c k (o,t) = p k (a) + tq k (a) , (k=l,2,3,4) , (A5) 

p.(») = i [• [st) b + ( "«+ y ) <,D + (~+ i)(ai>+bc)] , 

<h(“) = % [- (^)B • (^)«D • AD - BC] , 

Pj(«) = J [- <«-3)° ' + («+l)(AC-BD)] , 

q 2 («) = t t- a ‘ (^) A ’ (^)“ c - AC + BD J ’ 

p 3 (a) = 2 Pl (a) - ± , 
q 3 (a) = 2 qi (a) + , 

P 4 (a) = [(«-2)B - 3aD - 2AD - 2kBC + (/c+2)B(C 2 +D 2 ) + g (A 2 +B 2 )] , 

q 4 («) = [B + aD - B(C 2 +D 2 ) - g (A 2 +B 2 )] , (A6-13) 

Z = 2a[-B + /caD + 2 kAD + 2BC - B(C 2 +D 2 ) + (A 2 +B 2 )] , (A14) 

where A, B, C and D are functions of a and are defined by 
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A + iB = n 1 (a) , C + iD = m 1 (a) 


(A15,16) 


and nj and m 1 are given by (14a) and (12). 
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APPENDIX B 


The asymptotic expansions needed to examine the singular behavior 

of the kernel K(x,t) 


that 


For a — ►oo from Appendix A and equations (12), (14), (16) and (28) it can be 


R(a) = 2a + 4^“) a + °( a ' 3 ) * 


Q _z_ 1 

9 ~ 2S 


_3i k \i . ( 15 ~*)^ _L + 0(a' 5 )] 

1 +k [q 12(1 +/c) a 3 + ^ J 


“ a ‘ I + 2(fc) “ + 0( ^> ’ 


B - -2N 


1+k & ' 2 (fc) ^ + ° ( ^ )] ’ 


. , (2-«)/3 1 (/c 3 -2k 2 -3/c+2)/? 2 1 , (1+«)(2-k )(3-/c)^ _i_ n , 

= -1 - — y- a 4(TR ? + 8 a 3 + °V' ’ 


D 


- 1 Syl + (l+*)(2-/c)/? 2 X + (/c 4 -« 3 -5/c 2 -/c+4 )T ^ 0( JL)i 

” 2 \| 1+k [kP a + 2 a 2 + 4(l+/c) a 3 + ’ 


Z - -*1® (1-«)(1+«)W» - + 0(4)] , 


h l( a ) “ ^ + 0 (^) * 


a 


h 2 (a) = O(J^) , 
a 


h r.'i - 0 1 0 --*)P 2 1 , q/ 1 \ 

h 3<“) - - (ThT) + 5 + <V ’ 


shown 

(Bl) 

(B2) 

(B3) 

(B4) 

(B5) 

(B6) 

(B7) 

(B8) 

(B9) 

(BIO) 
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